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1 Introduction 

Statistical mechanical lattice models are called solvable if their associated 
Boltzmann weights satisfy the factorization or star-triangle equations of Mc- 
Guire Q, Yang [^| and Baxter ||. For such models the free energy per site 
and the one-point correlations in the thermodynamic limit are expressible in 
closed form [||. There exists a deep mathematical structure j|, || underly- 
ing these solvable models; and near critical points, a wider applicability than 
one would initially expect. This last phenomenon, called universality, has its 
mathematical roots in the strong law of large numbers and the central limit 
theorems of probability theory and its physical origins in critical phenonmena 
and conformal held theory. 
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The exact computation of n-point correlation functions is generally an open 
problem for most "solvable" models. In the special case of the 2D Ising model, 
the n-point functions (in the scaling limit) are expressible in terms of solutions 
to integrable differential equations || g] . The Wigner-Dyson theory of ran- 
dom matrices [|[ 10 is a second class of statistical models where integrable 



differential equations and n-point correlations (or more precisely, level-spacing 
distributions) are related. This paper reviews some of these relationships. 

In section 2 we define the basic objects of random matrix theory (RMT). 
In section 3 we recall the bulk scaling limit and the edge scaling limit and 
express various distribution functions in terms of Painleve transcendents. In 
section 4 we discuss the universality of these distribution functions. 



2 Random Matrix Models 

In the Gaussian models || [l0| , the probability density that the eigenvalues lie 
in infinitesimal intervals about the points X\, . . . , Xx is given by 

P N p(x 1) ...,x N ) = C N0 e-^ x l H\ Xj -x k f, (1) 

j<k 

where Cn/3 is a normalization constant and 

f 1 for GOE, 
(3 := I 2 for GUE, 
[ 4 for GSE. 

We recall that for (3 = 1 the matrices are N x N real symmetric, for (3 = 2 the 
matrices are NxN complex Hermitian, and for (3 = 4 the matrices are 2N x 2 N 
self-dual Hermitian matrices. (For (3 = A each eigenvalue has multiplicity two.) 
In each case the eigenvalues are real. 
In RMT the probabilities 

Enp^J) := J ■ ■ ■ J Pnp(xi, . . . , xn) dx\ ■ ■ ■ dxN, 

= probability no eigenvalues lie in J, 

are particularly interesting (To). The simplest choices of J are (a, b) and (f, oo). 
In the first instance the mixed second partial derivative E^p(0; J) with respect 
to a and b gives the spacing distribution between consecutive eigenvalues; and 
in the second case, FNp(t) := Epfp(0; (t, oo)) is the distribution function for 
the largest eigenvalue. 
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The important fact |ll], [T^, |l3j (sec |Q for simplified proofs) is that 
Enp(0',J) (or its square for (5 = 1 or 4) is expressible as a Fredholm de- 
terminant. For (3 — 2 the kernel is a scalar kernel acting on J whereas for 
(3 = 1, 4, Enp(0; J) 2 equals a Fredholm determinant of 2 x 2 matrix kernel 
acting on J. 



3 Painleve Representations 

The representation of E^p(0; J), or its square, as a Fredholm determinant 
of an operator with kernel of a special form permits the determinant to be 
expressed in terms of Painleve functions Jl5|, [l7], [l^, |l^] . 



3.1 Bulk Scaling Limit 

Denote the density of eigenvalues at the point xq by p(xo). It is customary 
in the limit N — > oo to scale distances so that the resulting density is one. 
Precisely, we define £ = pn{xq)(x — xq), xq independent of N, and consider the 
limit N — > oo, x — > scq, such that £ is fixed. By requiring xo to be independent 
of N, we are choosing a point in the "bulk" of the spectrum and are examining 
the local statistics of the eigenvalues in some small neighborhood of the point 
Xq. In this limit, and for j3 = 2, we are led to the Fredholm determinant of the 
operator on L 2 (0, s) whose kernel is the famous sine- kernel [O |o) 



Observe that the kernel is translationally invariant and independent of the 
point Xq. 

It is a result of Jimbo et al. 1 15 that 

dct (I - \K) = exp (^£ S dx^j , 

where a satisfies the differential equation 

(xa"f + 4 (xa' - a) (xa' -a+ {a') 2 ) = 

with the boundary condition a(x; A) ~ — > 0. (Here K is the operator 
whose kernel is the sine-kernel acting on L 2 (0,s) : ) The differential equation 
is the "a representation" of the Pv equation [jl5[ go) . For other proofs of this 
result see [El El 53. 
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If Ep(0; s) denotes the limiting value of En^(0; (— i, t)) in the bulk scaling 
limit with the scaled length of J set equal to s, then jlO|, |l6|, 



£i(0;s) = det(I-K+), 
E 2 (0;s) = det{I-K), 

£ 4 (0;s/2) = ^ (det (I - K + ) + det (I - K-)) 

where K± are the operators with kernels K(x,y)±K{— x, y). The det(7 — K±) 
can be expressed in terms of det(J — K) 15, |16| : 



det (7 - K±f = det (I - K) exp 



d 2 
dx 2 



log det (I - K) dx 



These formulas for Ep(0; s) are well adapted for producing graphs of Ep(0; s) 
and the level-spacing densities pp(s) = d 2 Eg(0;s)/ds 2 once one knows the 
large s asymptotics of a(s; 1) and Ep(0; s) (TGI, |25j| . 



3.2 Edge Scaling Limit 

The limiting law, called the Wigner semi-circle law, is well known 
1 



lim 

AT,,— oo 2(7 J N/3 



Pn (^(J^/Wb 



^VT^ \x\<i, 

\x\ > 1. 



Here a, <j/V2, a J \/2 (for ft = 1,2, 4, respectively) is the standard deviation of 
the Gaussian distribution in the off-diagonal elements and 

N,ft=l, 
N l3 ={ N, ft = 2, 

2N + 1, ft = 4. 

For the normalization here, a = 1/V2. Less known is that the distribution 
function of the largest eigenvalue satisfies [Efl : 



F 



N/3 



(2cr A 



N fj + x 



ifx<0, 

1 ifx>0, 



as N — > oo. The edge scaling variable, s, defined by 

a s 



t = 2a J N, 



w 



Table 1: The mean (up), standard deviation (erg), skewness (Sp) and kurtosis 
(Kp)oiFp. 





M/3 


op 


Sp 


Kp 


i 


-1.20653 


1.2680 


0.293 


0.165 


2 


-1.77109 


0.9018 


0.224 


0.093 


4 


-2.30688 


0.7195 


0.166 


0.050 



gives the scale of the fluctuations at the edge of the spectrum |g7|, £8|, [lTj . As 
N -> oo with s fixed @ [l8| 



Pv 2 (i) -» exp ^- jf 



(x - s)<z(x) 2 =:F 2 (s) (2) 



where q is the solution to the P// equation 

q" = s q + 2 g 3 

satisfying the condition 

g(s) ~ Ai(s) as oo, 

with Ai the Airy function. 

For (3 = 1,4 the results are 

Pjvi (*) 2 -> Pi(s) 2 = F 2 (s) exp ^- J q(x) dx 

and 

FAr 4 (i) 2 -» F 4 (s/V2) 2 = F 2 (s) cosh 2 Q jf q(x) 

Table 1 gives some statistics of Fp and Figure 1 shows the densities fp(s) — 
dFp/ds. 

4 Universality 

The universality of the distribution functions of RMT largely accounts for its 
success in a wide range of applications, see, e.g. |2?|. However, this universality 
has only been proved in a small number of cases. The evidence for universality 
is largely numerical. 
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Figure 1: The probability density fp of the largest eigenvalue, (3 — 1, 2,4. 
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There are two types of universality results. The first, and the easier to es- 
tablish rigorously, consists in modifying the random matrix model itself by, say, 
replacing the Gaussian x 2 /2 appearing in (|l]) by an "arbitrary" potential V(x). 
Changing the potential V changes the density of eigenvalues p, but several au- 
thors have established that the bulk scaling limit results in the sine-kernel and 
at soft edges one generically obtains the Airy universality class. (Fine tuning 
at the edge can result in different universality classes, see, e.g. |27|.) 

The second type of universality, and the one first envisioned by Wigner, 
asserts that for a classical, "fully" chaotic Hamiltonian the corresponding 
quantum system has a level spacing distribution equal to Pj3{s) in the bulk. 
(The symmetry class determines which ensemble.) A nice numerical example 
of this quantum chaos is Robnik's work ]3l| on chaotic billiards. 

Here we discuss three other examples of RMT universality. 

4.1 Zeros of the Riemann Zeta Function 

Work by Montgomery |?2| followed by extensive numerical calculations by 
Odlyzko ]33| on zeros of the Riemann zeta function have given convincing 
numerical evidence that the normalized consecutive spacings follow the GUE 
distribution, see Figure 2 (the GUE Hypothesis). Rudnick and Sarnak |34| 
have proved a restricted form of this hypothesis. 

4.2 Eigenvalues of Adjacency Matrices of Quasiperiodic 
Tilings 

The discovery of quasicrystals has made the study of statistical mechanical 
models whose underlying lattice is quasiperiodic of considerable interest to 
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Figure 2: Data for nearest neighbor spacings among 1,041,600 zeros of the 
Riemann zeta function near the 2 x 10 20 -th zero are plotted together with the 
GUE spacing density. Courtesy of Andrew Odlyzko p^ j. 

Nearest neighbor spacings among 
1 ,041 ,600 zeros near the 2*1 A 20-th zero 
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normalized spacing 



physicists. In particular, in order to understand transport properties, tight 
binding models have been defined on various quasiperiodic lattices. One such 
study by Zhong et al. f35|| defined a simplified tight binding model for the 
octagonal tiling of Ammann and Beenker. This quasiperiodic tiling consists 
of squares and rhombi with all edges of equal lengths (see, e.g., |3q] ) and has 
a Z?g symmetry around the central vertex. On this tiling the authors take 
as their Hamiltonian the adjacency matrix for the graph with free boundary 
conditions. The largest lattice they consider has 157,369 vertices. The matrix 
splits into ten blocks according to the irreducible representations of the dihedral 
group Dg. For each of these ten independent subspectra, they compare the 
empirical cumulative distribution of the normalized spacings of the consecutive 
eigenvalues with the integrated density 1% (s) — J p\ (x) dx where p\ is the 
GOE level spacing density. In Figure 3 we have reproduced a small portion of 
their data for one such subspectrum together with I\ . 

4.3 Distribution of the Length of the Longest Increasing 
Subsequence of Random Permutations 

An old problem going back to Ulam asks for the limiting behavior of the length 
of the longest increasing subsequence of a random permutation. Precisely, if 7r 
is a permutation of {1, 2, . . . , N}, we say that • ■ ■ , Tr(ik) is an increasing 
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Figure 3: Data for nearest neighbor, cumulative, normalized spacings of eigen- 
values of the adjacency matrix for a quasiperiodic octagonal tiling are plotted 
together with the GOE integrated spacing density ii(s). Data are from one 
independent subspectrum of a Ds-symmetric octagonal patch of a tiling with 
157,369 vertices. Courtesy of Zhong et al. |35[ ]. 
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subsequence in ir if ii < • • • < ik and n(ii) < • • • < n(ik)- Let £n{^) be the 
length of the longest increasing subsequence of n. We take each permutation 
to be equally likely thus making £n a random variable. The problem then is 
to understand the distribution of £n- This problem has a long history (see, 
e.g. [[36]]) but its connection with RMT is recent |37, In a recent paper, 
Baik, Deift and Johansson |36| have proved the following remarkable result: 



Theorem: Let Sn be the group of all permutations of N numbers with uni- 
form distribution and let £n{tt) be the length of the longest increasing subse- 
quence of 7T £ Sn- Let x be a random variable whose distribution function is 
F2, the distribution function for the largest eigenvalue in the GUE in the edge 
scaling limit (0). Then, as N — > 00, 



XN 



in distribution, i.e. 



lim Prob 



Xn 



iW6 



< 8 = F 2 (s) 



for all s £ R. 
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